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Abstract. In the present article we show that the energy spectrum of the one- dimensional Dirac equation, 
in the presence of an attractive vectorial delta potential, exhibits a resonant behavior when one includes an 
asymptotically spatially vanishing weak electric field associated with a hyperbolic tangent potential. We 
solve the Dirac equation in terms of Gauss hyper-geometric functions and show explicitly how the resonant 
behavior depends on the strength of the electric field evaluated at the support of the point interaction. We 
derive an approximate expression for the value of the resonances and compare the results calculated for 
the hyperbolic potential with those obtained for a linear perturbative potential. Finally, we characterize 
the resonances with the help of the phase shift and the Wigner delay time. 

PACS. 03.65.Pm Relativistic wave equations - 03.65.Ge Solutions of wave equations: bound states 



1 Introduction 

Spontaneous particle production in the presence of strong 
electromagnetic fields is undoubtedly one of the most in- 
teresting phenomena associated with the quantum vac- 
uum pQ. The study quantum effects associated with rela- 
tivistic particles in strong fields dates back to the pioneer- 
ing work of Klein [2] where he studied the reflection and 
. | penetration of electron waves on a potential barrier, ob- 
. taining the result that, for very strong potentials a large 
' number of electrons with negative energy penetrate into 
?3 \ the wall. This effect is called the Klein paradox. Sauter [5] 
found essentially the same results in the most general case 
when the barrier is smoothed. Particle creation by strong 
infinitely extended constant electric field, which is closely 
related to the Klein paradox, was predicted by Heisenberg 
and Euler 0). 

The study of supercritical effects and resonant par- 
ticle production by strong Coulomb-like potentials dates 
back to the pioneering works of Pieper and Greiner [5] 
and Gershtein and Zeldovich [5] where it was shown that 
spontaneous positron production was possible when two 
heavy bare nuclei with total charge larger than some crit- 
ical value Z c collided with each other. The critical Z c = 
1/a » 137 is the value for which the IS state of the hy- 
drogenlike atom with potential V = — Ze/r has energy 
E = —mc 2 . Supercriticality effects are based on sponta- 
neous positron emission induced by the presence of very 
strong attractive electric potentials. The energy level of 
an unoccupied bound state sinks into the negative energy 
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continuum. An electron of the Dirac sea is trapped by 
the potential, leaving a positron that escapes to infinity. 
The electric field responsible for supercritical effects has 
a strength of E ~ m 2 c 3 /(eh) and, in the region where 
the field acts, it makes a work larger than 2m e c 2 which 
is the value of the gap between the negative and positive 
energy continua. Such strong electric fields could be pro- 
duced in super-heavy nuclei, heavy-ion collisions [TJG] and 
in astrophysical phenomena. 

In order to get a deeper understanding of the mecha- 
nism responsible for resonant peaks appearing in the en- 
ergy spectrum when supercritical fields are present, and 
the role played by perturbative fields in the shape of res- 
onances, we proceed to work with a vector point inter- 
action in the presence of a homogeneous, asymptotically 
vanishing electric field [8,9,10 Point interactions poten- 
tials permit us to tackle, in a simple way, more complex 
short-ranged potentials. Among the advantages of work- 
ing with confining delta vector potentials we have the fact 
that, they only possess a single bound state and the treat- 
ment of the interaction reduces to a boundary condition. 
Bound states of the relativistic wave equation in the pres- 
ence of point interactions have been carefully discussed 
in the literature [Ill[T^l[T3"l[PT] . The one-dimensional Dirac 
equation in the presence of a vector point delta interaction 
has also been a subject of study in the search for supercrit- 
ical effects induced by attractive potentials. In this case, 
we see that a vectorial delta potential is strong enough to 
pull the bound state into the negative energy continuum 
E = —mc 2 P^llPollTfj] . nevertheless this supercritical state 
does not evolve to a real resonant state. Since we are inter- 
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cstcd in studying the mechanism of positron production 
by supercritical fields, we proceed to analyze the resonant 
behavior of the energy when a bound state dives into the 
negative continuum. This resonant behavior is associated 
with the appearance of simple poles of the resolvent on 
the unphysical sheet at a position very near the real axis 
|17j , which can be identified as positron states with a short 
mean life. 

In the present article, applying the idea developed by 
Titchmarsch [18] and Barut [19], we compute the energy 
spectrum of the one-dimensional Dirac equation in the 
presence of a vector Dirac delta interaction and a weak 
electric field associated with a linear and an hyperbolic po- 
tential. We find that, in both cases the energy spectrum 
exhibits a resonant behavior. This problem can be con- 
sidered the relativistic extension of the one-dimensional 
Stark effect [2T]l22ll23] . 

The paper is structured as follows: In section 2, we 
solve the one-dimensional Dirac equation in the presence 
of an attractive S potential and a hyperbolic tangent vec- 
tor potential. In section 3, we compute the energy reso- 
nances and show how they depend on the electric field 
strength evaluated at the delta interaction point. We also 
derive an approximate analytic expression for the energy 
resonances. In section 4, we characterize the resonance by 
computing the phase shift <j) and the Wigner time delay 
for the hyperbolic potential. Finally, in section 5 we sum- 
marize our conclusions. 



2 The one-dimensional Dirac equation 

The energy spectrum of the Dirac equation in the present 
of a delta potential has been studied by different authors 
in the literature. The attractive vectorial delta interaction 
is only is able to support a bound state which becomes 
a zero-momentum resonance for E = —m but it is not 
strong enough to sink the bound state into the Dirac see 
to create a time-decaying resonant state. In this section 
we proceed to solve the 1+1 Dirac equation in the pres- 
ence of an attractive vector point delta interaction and a 
potential V(x) which we will consider perturbative. The 
Dirac equation, expressed in natural units (U = c = 1) 
takes the form l20l 



%Y 



ieA, 



m # = 0, 



(1) 



where is the vector potential that in our case takes the 
form eA^ — (—gS(x) + V(x))6q, e is the charge, and m 
is the mass of the electron. The Dirac matrices 7 M satisfy 
the commutation relation {7^,7^} = 2r) fJ - lJ with rf v = 
diag(l, —1). Since we are working in 1+1 dimensions, we 
choose to work in a two-dimensional representation of the 
Dirac matrices 



7 



^3,7 



-l(T 2 . 



(2) 



Substituting the matrix representation ([2]) into equation 
([T]), and taking into account that the potential interaction 
does not depend on time, we obtain 



{-ivij- + (V(x) - E) + ma 3 }X(a;) = 0, (3) 



with & = 03X, and 



X(x) 



X 2 



(4) 



where E is the energy eigenvalue. The components of Eq. 
(UJ) satisfy the boundary conditions: 

Xi(0 + ) = Xi(O-) cos g- iX 2 (0~) sing, (5) 
X 2 (0 + ) = -iXi(O-) sing + X 2 (0-) cos g. 

The relations given by Eq. ([5]) describe an attractive delta 
vector potential interaction of strength g Q3] . 

Taking into account the gamma matrices representa- 
tion ([2]) , we obtain the system of equations 



(m + V(x) -E)Xx-i 
.dXx 



dX 2 
dx 



0, 



i— — + (m - V(x) +E)X 2 = 0, 
ax 

Introducing the functions fl\ and 2 

Xx = Ox + i0 2 , X 2 = Ox - i0 2 , 



(6) 
(7) 

(8) 

we see that the system of equations given by Eqs. ©-© 
reduces to 

^P- + i (V(x) -E)Ox- m0 2 = 0, (9) 
dx 

^ - i (V(x) -E)0 2 - raQ 1 = 0, (10) 
dx 

which is more tractable in the search of exact solutions. 
Using Eq. §§§ and Eq. (JTUJ) we obtain that ± (x) satisfies 
the second order differential equation 

^O.ix) + {*^^ + (V(x) - Ef - m 2 } O^x) = (11) 



2.1 Delta interaction and linear potential Ax 

We are interested in solving the Dirac equation |3]) in the 
presence of an attractive delta and a linear perturbative 
potential V(x) = A x, which corresponds to a constant 
electric field eE x = — A. The perturbation of the energy 
spectrum of a non-relativistic particle by a constant elec- 
tric field has been extensively discussed in the literature, 
mainly because to its relation to the Stark effect and the 
presence of energy resonances [17) . 

The solution to Eq. ([TTjl . for A < 0, satisfying the 
boundary condition of a growing oscillatory solution for 
x < and a damping solution for x > can be expressed 
in terms of the parabolic cylinder functions D v (x) [24] as 



0+(x) = AD 



0;{x)=BD p 




-ef [Xx — E) 



x > 0. 



ef (Xx- E) , x < 0. (12) 
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where D p are the parabolic cylinder functions [24 , p = 
A and B are constants to be determined using the 
boundary condition given by Eq. ([5]). The solution Q+(x) 
belongs to £ 2 (0, oo) with ImE > 0. The solution fi^(x) 
behaves asymptotically as an outgoing wave, a condition 
that defines a Siegert state [25] 

Inserting Eq. (fT2| into Eq. Q and using the recur- 
rence relations for the parabolic cylinder functions [24] , 
we obtain 



Q+{x)=A^=^D p ^\d-e^{Xx-E)\ (13) 



2A 

-%-k/a ( n \ 

Q-(x)=B y= D p ^i-Jje i! f(Xx-E)Uu) 



2X 



Using the boundary condition (|5|) we obtain that the en- 
ergy resonances satisfy the equation 

(A7(0)X+(0) - X 2 -(0)X+(0)) cos( 3 ) 
+i(X-(0)X+(0) - X-(0)X+(0)) sinQ?) = 0, (15) 

and with the help of the definition of X\ , and X 2 in terms 
of fix and fl 2 © , we obtain the result that the eigenvalue 
equation (| 1 5|) can be written as: 



(0)^(0) - e- 2i »/?r(0)fl 2 + (0) = (16) 

which, after substituting the expressions for fl\ and fl 2 in 
terms of the parabolic cylinder functions reduces to 



D a - 



t£ D 



-te 



-2 lgD 




0(17) 



The eigenvalue equation (|17|) permits us to obtain the en- 
ergy solutions to the Dirac equation Q in the presence 
of a vectorial delta potential and a linear potential Xx. 
For small values of A the linear potential plays the role 
of perturbing the Dirac hamiltonian with a vectorial delta 
interaction —gS(x) creating resonant states that are super- 
critical when g approaches ir and E — > —to, The solutions 
to Eq. (JTrJ) for small values of A should be obtained numer- 
ically. Fig. [1] shows QE versus 3tE for varying A showing 
a resonant behavior analogous to the one expected for the 
supercritical hydrogen atom. 



2.2 Delta interaction and a hyperbolic potential 

Now we proceed to solve the 1+1 Dirac equation in the 
presence of an attractive vector point interaction potential 
represented by eV(x) — —gS(x), and an asymptotically 
vanishing electric field associated with the potential 



V(x) = Atanh(foc). 



(18) 



The potential (|T8|) is asymptotically constant for large val- 
ues of x and has the linear potential as a limit for small 
values of k. 



E 
lD 
E 



0.05 




-1.02 

Re E/m 



Fig. 1. Behavior of Ss(E) against ^R(E) for the resonant energy 
states of a delta interaction with g — — 0.99-7T perturbed by 
a linear potential Aa; with < A < 0.03. 3f(I5) and -U(E) 
increase as the parameter A increases 

After substituting the potential (|18[) into Eq. (fTTj) . we 
obtain the second order differential equation 



d 2 f2i(x) 
dx 2 



In order to be able to identify the solutions of (|19|) ex- 
hibiting a damping asymptotic behavior as x — > +oo, we 



-{ikAsech 2 (kx)+(At&nh(kx)-E) 2 -m 2 }ni(x) = 

(19) 



introduce, for x > 0, the new variable z 



-Ikx 



We see 



that fi\(z) satisfies the second order differential equation 
,d 2 I2i(z) 



Ak 2 z 2 



ik 2 z dQl ^ 
dz 2 dz 



■ AzAk 

'WW 



Q 1 {z)=Q. (20) 



For k > and A < 0, the solution of (f2"0"l) . vanishing as 
x — ► oo , with dlE < —to and QE > 0, can be expressed 
in terms of Gauss hyper-geometric functions 2 Fi (a, 6, c, z) 
\M as 



n+(z) =Cx{z- ly^z^F, (c, d, /; z) , (21) 



where 



Jm 2 -{E- A) 2 Jm 2 - (E + A) 2 
a= 2k ' b= 2k ' 

c = a — 6 . d = a + b(E, A 7 m 7 k) — i — , 

k k 

f = 1 + 2a. (22) 

Using Eq. © we can obtain the expression for (z) In 
order to solve Eq. (fl9|) for a: < 0, we introduce the new 



variable z 
equation 



,.,2kx 



obtaining the following differential 



ik^z 



+ 



■ lr . 2 d 2 f2i(z) 



dl 2 

; AzAk 

'WW 



dz 

\ 2 

- E) -to 2 



Qx{z) =0 (23) 
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whose solution, exhibiting the behavior of an increasing 
reflected wave as x — ► — do. with 'RE < —m and 'HE > 0, 
is 



fir (2) = Ca(2f- 2F1 (ff, h,p;2) , 



(24) 



with 



3 = -6 
p= -1- 26 



a — , h = —b — a — , 



(25) 



The energy eigenvalues equation corresponding to the Dirac 
equation Eq. fl} , with the boundary conditions {JSJ can be 
obtained after substituting fl\ and J?2 into Eq. (fT6|) . When 
the potential V(x) vanishes at x — 0, such as in the linear 
potential or hyperbolic cases, the eigenvalue equation (|16l) 
can be written in terms of Q~ (x) and (x) as 



dD ^ x) nt{x)-e-^ dQ l {x) Q;{x) 
-iE(l - e- 2tg )f2-(x)f2+(x) = 



(26) 



3 Approximate solutions 



It is not straightforward to obtain an approximate expres- 
sion for the energy eigenvalues of Eq. (|17|) in the vicinity 
of E = — m because the expansion parameter A appears 
in Eq. p7[) in the argument and order of the parabolic 
cylinder functions. It is also not possible to try to apply 
perturbation theory to find the complex energy eigenval- 
ues. In this section, we derive an approximate solution to 
Eq.Q in terms of Airy functions for the linear potential 
A x and for the hyperbolic potential A tanh(fca;) for small 
A and k. 

Since we are interested in studying Eq. (jTTJ) for very 
small values of A and search solutions of Eq. (fTS)) with 
respect to the resonant energy value E = — m, for the 
linear potential V(x) = Ax, we approximate in Eq. (jTTJl , 
(A x — E) 2 by —2XEx + E 2 , obtaining in this way the 
approximate differential equation 



QAx) + \iX-2XEx 

dx 2 v ; 1 



E 2 -to 2 }Q 1 {x) = 0. (27) 



We proceed to solve Eq. (|2T|) demanding the solutions to 
satisfy the resonance asymptotic conditions, that is, we 
choose damping solutions for x > and diverging oscil- 
lating functions for x < 0. Those solutions are: 



n+{x) = AAi 



Q-{x) =BCi A 



iX + E 2 - m 2 - 2EX: 



(2EX) 2/3 
IX + E 2 - to 2 - 



2EX: 



(28) 



(29) 



{2EX) 2/3 

where A and B are constants and Ai(x) and Ci + {x) = 
Bi(x) + iAi(x) are the Airy functions [24] . 

Using Eq. (|2"8")l and Eq. (j2"9")l and the eigenvalue equa- 
tion (|26p . we obtain that the approximate spectral equa- 
tion takes the form 



{2EXy /d ze 



- + 



(30) 



2i Ai (a) (i{2EX) 1/a Ci+' {-) + ECi+ (-)) sin g = 



—E 



Where - " (2 EX) 2 '- ■ 

For the hyperbolic potential yltan(fca;) we also have 
that since the parameters A and k are small, we can ap- 
proximate the expression (Atanh(kx)~E) 2 by ~2AkxE+ 
E 2 in Eq. (|19p , obtaining in this way the approximate dif- 
ferential equation for J?j , 



dx 



-n i {x) + {iAk-2AkEx + E 2 - to 2 } n x (x) = 0(31) 



which , after making the identification Ak — > A reduces to 
Eq. (|27[) . therefore the approximate eigenvalue equation is 
given by Eq. ([3DJ) with A = Ak. 

In order to show the accuracy of the approximations 
(|3"Tj) and P?| , we can compare the solutions l[28p and ([2^)1 
with the exact solutions given by Eq. (p~2|) . Fig. 2 shows 
that, for small values of A, the solutions (|28|) and (|29|) are 
a good approximation to the parabolic cylinder functions 
([TZ]), and the Gauss functions ((2TJ) and (J24|) for small A 
and k. A comparison between the resonant energy values 
obtained after solving Eq. (|30|) and Eq. (|17p. shows that, 
for A = —0.01, the eigenvalue equation ([30]) gives resonant 
energies with an error smaller than 0.1% . For A = —0.01 
the solution to Eq. fTfl gives E = -1.016964 + 0.049316/ 
and the approximate solution obtained after solving Eq. 
(J3DJ) is E = -1.016139 + 0.0499657. Eq. ^ gives more 
accurate energy eigenvalues as A — > 0. The advantage of 
using Eq. (|30|) for calculating the energy eigenvalues for 
small values of A, instead of solving Eq. (jTTJ) , lies in the fact 
that, for large values of the index p, the parabolic cylinder 
functions exhibit a divergent behavior that makes it trou- 
blesome to obtain accurate values for the energy spectrum. 
For values of |A| < 0.01 we made use of the approximate 
solutions in order calculate the resonance energies. 




Fig. 2. Comparison between the exact solutions for Hi(x) 
(dashed line) and the approximate expressions (solid line) for a 
particle in an attractive delta perturbed by a linear potential 
Aa; with A = —0.1. Left and right figures depict respectively 
on the left and right of the origin x = 
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The approximation (|27|) to Eq. (fl9|) is valid when A 
and k are small with respect to m. For large values of 
k the hyperbolic potential approaches a step potential, 
which is not able to sink the delta bound energy state 
negative energy continuum. Fig. [3] shows the dependence 
of the resonant energy as the slope parameter k of the 
hyperbolic potential increases. 



°- 2 5| 




Re E/m 



Fig. 3. Behavior of the resonant energy for a delta attrac- 
tive interaction with g — — 0.997T perturbed by the potential 
Atanh(kx) with A = — O.f and < k < 0.5. For small values 
of k E is close to — f and Q{E) increases as k increases. 



4 Phase shifts. 

Since the perturbative potential yltanh(fca;) is asymptoti- 
cally constant as x — > ±oo, and, the perturbative electric 
field Aksech(kx) 2 / e vanishes for large values of x, we can 
use of the scattering formalism in order to compute the 
phase shift 0, magnitude that will help us compute the in- 
fluence of the perturbative hyperbolic interaction on the 
energy spectrum of the delta interaction. Analogous to 
the non-relativistic case [2T] , the particle is constrained to 
move in a half line and the reflection amplitude r satisfies 
|r| = 1. The phase shift <j) corresponds to the argument of 
r as <f> = — arg(r). 

With the help of the boundary condition associated 
with the vectorial delta interaction —g8{x), we proceed to 
calculate the reflection amplitude r 

tX{ = (X\ + rX{) cos g - {X\ + rX T 2 ) sing (32) 

tX\ = -i {X{ + rX{) sing + (X\ + rX r 2 ) cos.g, (33) 

where t and X 1 are respectively the transmission coeffi- 
cient and the transmitted wave, X % denotes the incident 
wave. From Eq. (|32|) and Eq. (f33|) . we see that 

T ~ {X{Xl - X\X r 2 ) cosg + i (X{X{ - A^AT*) sin g ' 

(34) 



The solution to Eq. (| 19(1 , behaving as wave incoming from 
the left is 

n^ ln {z)=D 1 (z-l)-^z b 2 F 1 (q,t,l + 2b;T), (35) 
where 

q = -a + b- l 4- (36) 
k 

t = -a + b- 1 -^- (37) 
k 

The component Q 2 m (z) can be obtained using Eq. ©. 
Figure [4] shows the behavior of the phase shift </> against 
the real part of the energy E for k = O.f and A = —O.f. 
It can be observed that the phase shift increases and has 
a jump of 7r as it approaches the value tt/2 and E is close 
to the real value of the resonance E res = —1.01837 + 
0.05058/ This result indicates that we are in the presence 
of a bound state that dissolves itself into the negative 
energy continuum. 




Fig. 4. Behavior of the phase shift against the energy E for 
an attractive delta perturbed by the potential A tanh(fca:) with 
A = —O.f and k — O.f. The phase shift suffers a jump of n 
close to the real part of the resonant energy E rea . 

The slope of the curve shown in Figure 0] is related to 
the Wigner time delay r and can be estimated using the 
Heisenberg principle [26j 

7T 1 , „ 

r ^ W{ = ~ (38) 

\dE\E=Eo 

The Wigner time r permits us to estimate how long the 
bound state lives before it decays into the negative energy 
continuum. For A = —O.f and k — O.f we obtain that 
t = 0.04622, a value which is very close to ^s(E res ). 
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The inverse of the imaginary part of the energy E gives 
a measure of the mean of the-life time associated with the 
energy E [26l|27] , a state that can be interpreted as a 
positron, i.e a hole with a positive charge, created in the 
Dirac sea with a finite life-time. 



5 Discussion of the results 

The vector point interaction vectorial —gS(x) described 
by the boundary conditions ([5]) produces zero-momentum 
supercritical states for g = tt but they are not complex 
energy states associated with resonant states with a Breit- 
Wigner profile. 

In this article we have analyzed the energy spectrum 
of the one-dimensional Dirac equation in the presence 
of an attractive vectorial delta point interaction when 
one introduces perturbative potentials of the form Xx and 
A tanh(fcx) with small A and k, corresponding respectively 
to a constant electric field and an an asymptotically van- 
ishing electric field that reaches its maximal strength Ak 
at x = 0. The presence of the perturbative electric field in- 
duces the appearance of a resonant energy state in the one- 
dimensional Dirac equation. This phenomenon is analo- 
gous to the Stark effect [IT] . 

Since the vector interaction —gS(x) is not strong enough 
to produce a real resonance pT jlTMTB] . we have that the 
presence of the perturbative potential plays a crucial role 
in the appearance of the resonant energy state. We have 
studied and discussed the presence of supercriticality in 
two different ways: first, studying the resonant behavior 
of the energy spectrum in the vicinity of E = — m, and 
second, analyzing the behaviour of the phase shift as a 
function of the energy. 

With the help of the phase shift and the Wigner time- 
delay, we have shown that the mean life t of the resonant 
bound state decreases as the strength of the perturba- 
tive potential increases, a result which is analogous to the 
one observed when a perturbative constant electric field 
is considered. 

The perturbative potential Ax diverges as x — > ±oo. 
This asymptotic behavior prevents us from obtaining asymp- 
totically free solutions, therefore we cannot use standard 
scattering theory for obtaining the energy resonances [27] . 
In this sense, the exact and approximate approaches pre- 
sented in this article permits us to obtain the value of the 
resonances in cases where the potential does not go to zero 
at infinity. 
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